Greedy routing on networks of mobile agents 
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In this paper, we design a greedy routing on networks of mobile agents. In the greedy routing 
algorithm, every time step a packet in agent i is delivered to the agent j whose distance from the 
destination is shortest among searched neighbors of agent i. Based on the greedy routing, we study 
the traffic dynamics and traffic-driven epidemic spreading on networks of mobile agents. We find that 
the transportation capacity of networks and the epidemic threshold increase as the communication 
radius increases. For moderate moving speed, the transportation capacity of networks is the highest 
and the epidemic threshold maintains a large value. These results can help controUing the traffic 
congestion and epidemic spreading on mobile networks. 
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I. INTRODUCTION 

Due to the increasing importance of communication 
networks such as the Internet [l| and networks of mobile 
phone users Q in modern society, the traffic of informa- 
tion flows Q and the spreading of computer viruses [J] in 
these networks have attracted more and more attention. 

Traffic dynamics concern mainly how to effectively de- 
hver information packets from source to destination on a 
network and avoid traffic congestion. The study of epi- 
demic spreading focuses on the forecast and control of 
computer viruses. For a long time, the two types of dy- 
namical processes were investigated separately. Recently, 
Meloni et al. f5l| introduced a theoretical approach to in- 
corporate traffic dynamics in virus spreading. In partic- 
ular, they cast the susceptible-infected-susceptible (SIS) 
model [6[ in a flow scenario where contagion is carried 
by information packets traveling across the network. A 
susceptible node is more likely to be infected if it re- 
ceives more packets from infected neighbors. They found 
that the epidemic threshold depends on the traffic and 
decreases as traffic flow increases. 

In communication networks, information packets are 
forwarded from sources to destinations by specific routing 
protocol. To enhance the the transportation capacity of 
communication networks, researchers have designed vari- 
ous routing algorithms, including the shortest path_0,[lj, 
the integration of static and dynamic information [9| , the 
local routing [loj . the efficient routing [ll|, the greedy 
algorithm [12], and so on. It has been found that the 
routing algorithm also plays an important role in traffic- 
driven epidemic spreading [l^ [l^ . 

Previous studies about traffic dynamics usually focus 
on static networks, where nodes are motionless and links 
among nodes keep fixed. In our recent paper [isj . we have 
studied the transportation dynamics on networks of mo- 
bile agents. We assume that agents move on a plane and 
the searching area of a agent i is a circle centered at i . In 
that paper, information packets were delivered according 



to random routing algorithm, that is, a packet in agent 
i is forwarded to a randomly chosen agent j from agent 
j's searching area. The random routing can be applied in 
the case where moving agents cannot obtain the informa- 
tion about other agents' positions. However, if an agent 
can know other agents' positions, the random routing is 
not an effective algorithm for a packet to quickly reach 
the destination. Utilizing the information of agents' po- 
sitions, we now propose a greedy routing on networks of 
mobile agents. In the greedy routing, every time step 
a packet in agent i is delivered to the agent j whose 
distance from the destination is shortest among agent 
i's searched neighbors. We have found that the greedy 
routing markedly enhances the transportation capacity of 
networks compared with random routing. Based on the 
greedy routing, we study how moving speed and commu- 
nication radius affects traffic dynamics and traffic-driven 
epidemic spreading. 

The paper is organized as follows. In Sec. |lll we de- 
scribe our model in terms of the greedy routing and epi- 
demic spreading. Results on the traffic dynamics and the 
epidemic spreading are presented in Sec. Illll and Sec. IIVI 
respectively. A brief conclusion is given in Sec. El 



II. MODEL 

In our model, N agents move on a square-shaped cell of 
size L with periodic boundary conditions. Agents change 
their directions of motion 9 as time evolves, but the mov- 
ing speed V is the same for all agents. Initially, agents are 
randomly distributed on the cell. After each time step, 
the position and moving direction of an arbitrary agent 
i are updated according to 



Xi{t + 1) = Xi{t) + vcos9i{t), 



y.,{t + 1) = y,{t) + 
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where Xi{t) and yi{t) are the coordinates of the agent 
at time t, and ^'i is an iV-independent random variable 
uniformly distributed in the interval [— 7r,7r]. 

Each agent has the same communication radius r. Two 
agents can communicate with each other if the distance 
between them is less than r. We define an agent i's neigh- 
bors as agents who are within agent z's communication 
area. At each time step, there are R packets generated 
in the system, with randomly chosen source and destina- 
tion agents, and each agent can deliver at most C packet 
toward its destination. To transport a packet whose des- 
tination, an agent performs a local search within a circle 
of radius r. If the packet's destination is found within the 
searched area, it will be delivered directly to the destina- 
tion. Otherwise, the packet is forwarded to a according 
to the greedy routing, which is defined as follows. 

At time t, agent i has a packet whose destination is 
agent j. Agent k is one of z's neighbors at that moment. 
The distance between agent k and agent j is 



dMt) = ^J[xAt)-Xk{tW + Mt)^yk{tW. (4) 

Agent i will deliver the packet to the agent k whose dis- 
tance djk{t) is shortest among agent «'s neighbors. 

The queue length of each agent is assumed to be un- 
limited and the first-in-first-out principle holds for the 
queue. Once a packet reaches its destination, it will be 
removed from the system. 

After a transient time, the total number of delivered 
packets at each time will reach to a steady value, then 
an initial fraction of agents po is set to be infected (e.g., 
we set Po = 0.1 in numerical experiments). The infection 
spreads in the network through packet exchanges. For 
example, at time t agent i is infected and a packet is 
traveling from agent i to a susceptible agent j , then at the 
next time step, agent j will be infected with probability 
/3. The infected agents are recovered at rate n (we set 
/i = 1 in this paper). 



III. TRAFFIC DYNAMICS 

In this section, we study traffic dynamics based on the 
greedy routing. We set the number of agents N = 1500, 
the size of the square region L = 10 and the delivering 
ability of each agent C = 1 in this section. 

To characterize the transportation capacity of a net- 
work, we exploit the order parameter rj introduced in 
Ref. 0: 
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where AA'^n 



Np(t + At) - Np(t), (• • •) indicates the 
average over a time window of width At, and Np{t) rep- 
resents the total number of information packets in the 
whole network at time t. 
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FIG. 1: (Color online) The order parameter rj{R) as a function 
of the packet-generation rate R for random routing and greedy 
routing. For both routing algorithms, the moving speed v = 
0.1 and the communication radius r — 1. ri{R) is obtained by 
averaging over 5 x 10* time steps after disregarding 5 x 10'^ 
initial steps as transients. Each data point results from an 
average over 50 different realizations. 



As the packet-generation rate R is increased through 
a critical value of i?c, a transition occurs from free flow 
to congestion. For R < Rc, there is a balance between 
the number of generated and that of removed packets so 
that (ANp) = 0, leading to ri{R) = 0. In contrast, for 
R > Rc, congestion occurs and packets will accumulate 
in the system, resulting in a positive value of ri(R). The 
transportation capacity of a network can thus be charac- 
terized by the critical value Rc- 

Figure [1] shows that the order parameter ?y(i?) as a 
function of R for random routing and greedy routing. 
One can see that the critical value Rc for greedy routing 
is much larger than that of random routing when other 
parameters are the same. For the moving speed v — 0.1 
and the communication radius r — 1, the critical value 
Rc for random routing and greedy routing is about 8 and 
270 respectively. This result shows that the greedy rout- 
ing can greatly enhance the transportation capacity of 
networks in comparison with that of the random rout- 
ing. In the following, we will detailedly study the traffic 
dynamics based on the greedy routing. 

Figure shows that the order parameter ri{R) as a 
function of R for different values of v and r. One can 
see that both the moving speed and the communication 
radius affect the onset of traffic congestion. It is inter- 
esting to find that there are two different kinds of phase 
transition from free flow to congestion. When the mov- 
ing speed is very small (i.e., v — 0.001), the jamming 
transition is a second-order phase transition. When the 
moving speed is large enough (i.e., v = 0.1 and v = 1), 
the appearance of a congested phase belongs to a first- 
order phase transition. 
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FIG. 2: (Color online) (a) The order parameter rj{R) as a 
function of the packet-generation rate R for different values 
of the moving speed v. The communication radius r = 1. 
(b) The order parameter ri{R) vs R for different values of r. 
The moving speed v — 0.1. Each data point results from an 
average over 50 different realizations. 
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FIG. 3: (Color online) (a) The critical value Rc as a function 
of the moving speed v for different values of the communica- 
tion radius r. (b) The dependence of Rc on r for different 
values of v. Each data point results from an average over 50 
different realizations. 
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FIG. 4: (Color online) (a) The average traveling time (T) as a 
function of the packet-generation rate R (< Rc) for different 
values of the moving speed v. The communication radius 
r = \. (b) The dependence of (T) on R for different values of 
r. The moving speed v = 0.1. 



Since the transportation capacity of networks is char- 
acterized by the critical value i?c, it is of particular in- 
terest for us to study how i?c is affected by the moving 
speed V and the communication radius r. Figure 
shows the dependence of i?c on v for different values of r. 
One can observe a nonmonotonic behavior. For a fixed 
value of r, the largest value of i?c is obtained when the 
moving speed v is moderate (about 0.1). This result is 
different from that of Ref. [l^, where i?c increases as v 
increases for random routing. Figure |3l^b) shows the de- 
pendence of i?c on r for different values of v. One can 
see that, i?c increases as the communication radius r in- 
creases regardless of the values of v. 

To understand how the moving speed and communi- 
cation radius affect the transportation capacity of net- 
works, we study the average traveling time (T) and the 
load distribution P(n) in the next. 

The average traveling time (T) is defined as the av- 
erage time steps for a packet traveling from source to 
destination. We investigate the average traveling time 
(T) as a function of the packet-generation rate R in the 
free flow state (i? < R^. From Fig. 01 one can find that 
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FIG. 5: (Color online) The load distribution P(n) vs n for difTerent values of the moving speed v and the communication radius 
r. The packet-generation rate R = 100. From (a) to (d), v is 0.001, 0.01, 0.1, 0.8 respectively and r is fixed to be 1. From (e) 
to (h), r is 0.8, 1,1, 2, 3 respectively and v is fixed to be 0.01. 



(T) increases as R increases. Due to the randomness of 
movement and packet generation, the number of packets 
in an agent's queue fluctuates as time evolves. At a time, 
if the number of packets in an agent's queue excesses the 
agent's delivering ability, extra packets must wait for be- 
ing delivered in the next time step. As R increases, the 
above phenomenon of transient congestion occurs more 
frequently, leading to a longer (T). For a fixed value 
of R, we find that (T) keeps almost unchanged when the 
moving speed increases from 0.001 to 0.1 [see Fig.|31Ja)]. 
But as V continually increases, (T) increases. As is shown 
in Fig. lUa), for the same R, the average traveling time 
(T) increases as v increases from 0.1 to 0.7. Figure IDJb) 
shows the average traveling time (T) as function of R 
for different values of the communication radius r. One 
can see that, for the same value of R, (T) decreases as r 
increases. 

Figure [5] shows the load distribution P{n) vs n for dif- 
ferent values of v and r, where n represents the load and 
P{n) is the probability that an agent has the load n. We 
define an agent's load Ui as: 



AT ' 



(6) 



where riiit) is the number of packets staying in i's queue 
at time t, AT is period of time (we set AT = 5 x 10'') 
and the sum runs over a period of time AT. An agent's 
load Hi reflects the average number of packets staying in 
i's queue at a time. 



Figures [Sl^a)-(d) show the load distribution P{n) for 
different values of the moving speed v when the packet- 
generation rate R = 100 and the communication radius 
r = 1. One can see that, P(n) approximately displays the 
Poisson distribution when v is very small (i.e., v ~ 0.001 
and V — 0.01). As v increases (i.e., v — 0.1 and v — 0.8), 
the load distribution becomes highly homogeneous and 
all agents have almost the same load. The highest load 
of an agent decreases as v increases from v = 0.001 to 

V = 0.1 [see Figs. [Sl[a)-(c)]. It is noted that for the 
same packet-generation rate, the average traveling time 
(T) keeps almost unchanged for v = 0.001 and v = 0.1 
[see Fig. IH^a)]. Thus, the change of the highest load for 

V = 0.001 and z; = 0.1 is caused by the alteration of 
network structures. As the moving speed increases, each 
agent's neighbors change more frequently and the topol- 
ogy of network turns from quasistatic structure to dy- 
namical structure. When v increases from 0.1 to 0.8, the 
average traveling time (T) increases and the load of an 
agent enhances [see Figs. [5] (c) and (d)]. Since all agents 
have the same delivering ability C, the transportation 
capacity of the whole network is determined by the high- 
est load of an agent. The increase of the highest load 
indicates that the decrease of the network's transporta- 
tion capacity. Combining Fig. Ufa) with Figs. [Ha)-(d), 
we can understand why the highest transportation ca- 
pacity of networks is reached at the moderate value of 
the moving speed. 

Figures [He)-(h) show the load distribution P{n) for 
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FIG. 6: (Color online) (a) The average traveling time (T) as 
a function of the moving speed v for different values of the 
communication radius r. (b) The dependence of (T) on r for 
different values of v. 
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FIG. 7: (Color online) (a) Density of infected agents p as a 
function of the spreading rate /3 for different values of the mov- 
ing speed V. The communication radius r = 1.4. (b) Density 
of infected agents p as a function of jS for different values of 
r. The moving speed v = 0.1. The packet-generation rate is 
R — 4000 and each agent's delivering ability is infinite. Each 
curve results from an average over 50 different realizations. 



different values of r when v is fixed to be 0.01. One can 
observe that, P{n) changes from the approximate Pois- 
son distribution to highly homogeneous distribution as 
the communication radius r increases. Besides, the high- 
est load of an agent decreases as r increases. According 
to Fig.lUb), for the same packet-generation rate, the av- 
erage travehng time (T) decreases as r increases. The 
decrease of (T) shortens the time steps that a packet 
stays in the system and relieves the traffic load, leading 
to the enhancement of the network's transportation ca- 
pacity. Thus the phenomenon observed in Fig. |3^b) can 
be explained. 



IV. EPIDEMIC SPREADING 

In this section, we study how the greedy routing affects 
the traffic-driven epidemic spreading. We first consider 
the case where each agent's delivering ability is infinite, 
C — oo, so that traffic congestion will not occur in the 
network. 



In the case of infinite delivering ability, the average 
traveling time (T) is independent of packet-generation 
rate R. Figure [S] shows the dependence of (T) on the 
moving speed v and the communication radius r. For a 
fixed value of r, (T) keeps almost unchanged for v < 0.3 
but increases as v continually increases [see Fig. El^a)]. 
For a fixed value of v, (T) decreases as r increases [see 
Fig.EIb)]. 

Figure [7] shows the density of infected agents p as a 
function of the spreading rate /3 for different values of 
the moving speed v and the communication radius r. We 
observe that, there exists an epidemic threshold /3c, be- 
yond which the density of infected agents p is nonzero 
and increases as (3 is increased. For /3 < /3c, the epidemic 
dies out and p = 0. From Fig. [3 one can find that both 
V and r can affect the density of infected agents and the 
epidemic threshold. Large values of v and small values 
of r promote the epidemic spreading. 

Figure IHI^a) shows the epidemic threshold /3c as a func- 
tion of the moving speed v for different values of the 
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FIG. 8: (Color online) (a) Epidemic threshold /3c as a function 
of the moving speed v for different values of the communica- 
tion radius r. (b) Epidemic threshold Pc vs r for different 
values of v. The packet-generation rate R — 4000 and each 
agent's delivering ability is infinite. Each data point results 
from an average over 50 different realizations. The curves are 
theoretical predictions from Eq. ([9}. 



communication radius r. One can see that, for a fixed 
value of r, f3c keeps almost unchanged for d < 0.3 but 
decreases as v continually increases. Figure EUb) shows 
the dependence of /3c on r for different values of v. One 
can observe that /3c increases as r increases. In the fol- 
lowing, we will provide a theoretical analysis to calculate 
the epidemic threshold /3c. 

Since all agents have the same communication radius 
r and the load distribution P{n) follows the approximate 
Poisson distribution or highly homogeneous distribution, 
the mean-field theory can be applied. At each time step, 
the average number of packets that an agent delivers is 
R{T)/N. Thus the rate equation for the epidemic dy- 
namics can be written as 



dpjt) 
dt 



= -Pit) 



R{T) 
N 



/3p{t)[l-p(t)]. 



(7) 
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FIG. 9: (Color online) Epidemic threshold Pc as a function of 
the packet-generation rate R for C = 10 and C — > oo. The 
moving speed v = 0.5 and the communication radius r = 1. 
For C = 10, the critical packet-generating rate Rc is about 
2000. Each data point results from an average over 50 dif- 
ferent realizations. The solid and dash line is the theoretical 
prediction from Eq. (|9} and Eq. (fTT|) . respectively. 



we obtain 

p[-i + ^/?(i-p)] = o. 

From Eq.([8l), we get the epidemic threshold 



R{T)- 



(8) 



(9) 



After imposing the stationarity condition dp{t)/dt — 0, 



The comparison between numerical and theoretical val- 
ues of /3c is shown in Fig. |8] (we use the numerical values 
of (T) in Eq. © because of the difficulty in calculating 
(T) theoretically). From Fig. [51 we find that the theo- 
retical prediction of /3c is in good agreement with that of 
simulation result. 

Next, we consider the finite delivery ability. In 
this case, traffic congestion occurs when the packet- 
generation rate R is above a critical value Rc ■ Figure [9] 
shows the epidemic threshold /3c as a function of R for 
C = 10 and C — !■ cxd. One can see that, /3c scales in- 
versely with the i? for C CXD, as predicted by Eq. 
For C = 10, /3c decreases to a steady value as R increases. 
When R < Rc, f3c is almost the same for finite and in- 
finite delivery ability. However, for R > Rc, we observe 
that the value of /3c is larger for C = 10 than that of 
C oo, indicating that traffic congestion can suppress 
the spreading of disease. Similar phenomena were also 
been found in Rcfs. 5, 14]. 

In the case of finite delivery ability, the lower limit of 
the epidemic threshold can be estimated as follows. 

For sufficiently large values of R, all agents will become 
congested and each agent can deliver only C packets at 
each time step. Thus Eq. ^ must be revised as 



7 



dpjt) 
dt 



-p{t)+cppm-p{t)]. 



(10) 



Imposing the stationarity condition dp{t) /dt = and we 
obtain the epidemic threshold 



I3c = 



1 

C' 



(11) 



for infinite delivering ability, the epidemic threshold de- 
creases as the moving speed increases and increases as the 
communication radius increases. Fourthly, traffic conges- 
tion can suppress epidemic spreading. Since the study of 
networks of mobile agents has received increasing atten- 
tion in recent years, our results are valuable for under- 
standing the structure and dynamics of mobile networks. 



V. CONCLUSION 

In conclusion, we studied how greedy routing aflFects 
traffic dynamics and traffic-driven epidemic spreading on 
networks of mobile agents. Our main findings are the 
following. Firstly, for finite delivering ability, the trans- 
portation capacity of networks peaks at the moderate 
moving speed and increases as the communication ra- 
dius increases. Secondly, the average traveling time of a 
packet increases as the moving speed increases but de- 
creases as the communication radius increases. Thirdly, 
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